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ABSTRACT 

The recent discovery of the first known accretion-powered millisecond pulsar with 
the Rossi X-Ray Timing Explorer provides the first direct probe of the interaction of 
an accretion disk with the magnetic field of a weakly magnetic (B ^ 10 10 G) neutron 
star. We demonstrate that the presence of coherent pulsations from a weakly magnetic 
neutron star over a wide range of accretion rates places strong constraints on models 
of the disk-magnetosphere interaction. We argue that the simple M 3 / 7 scaling law for 
the Keplerian frequency at the magnetic interaction radius, widely used to model disk 
accretion onto magnetic stars, is not consistent with observations of SAX J1808. 4-3658 
for most proposed equations of state for stable neutron stars. We show that the usually 
neglected effects of multipole magnetic moments, radiation drag forces, and general 
relativity must be considered when modeling such weakly magnetic systems. 

Using only very general assumptions, we obtain a robust estimate of fj, ~ (1 — 10) x 
10 26 G cm 3 for the dipole magnetic moment of SAX J1808. 4-3658, implying a surface 
dipole field of ~ 10 8 — 10 9 G at the stellar equator. We therefore infer that after the 
end of its accretion phase, this source will become a normal millisecond radio pulsar. 
Finally, we compare the physical properties of this pulsar to those of the non-pulsing, 
weakly magnetic neutron stars in low-mass X-ray binaries and argue that the absence of 
coherent pulsations from the latter does not necessarily imply that these neutron stars 
have significantly different magnetic field strengths from SAX J1808. 4-3658. 



Subject headings: accretion, accretion disks — pulsars: individual: SAX J1808. 4-3658 
— stars: neutron — X-rays: stars 
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1. INTRODUCTION 

The basic framework for understanding accretion-powered pulsars emerged soon after their 
discovery (Giacconi et al. 1971). These systems are rotating neutron stars accreting matter from 
a binary companion, with magnetic fields strong enough to disrupt the accretion flow above the 
stellar surface (Pringle & Rees 1972; Davidson & Ostriker 1973; Lamb, Pethick, & Pines 1973). 
When threaded by the stellar magnetic field, the accreting gas is brought into corotation with the 
star and is channeled along field lines to the polar caps, releasing its potential and kinetic energy 
mostly in X-rays. The rotation of these hot spots through our line of sight produces X-ray pulses 
at the spin frequency of the neutron star. Most of the ss 50 accretion-powered pulsars have spin 
frequencies and inferred dipole magnetic field strengths in the range ~1 — 10 3 s and ~10 n — 10 13 G, 
respectively (White, Nagase, & Parmar 1995). 

In strongly magnetic (B ^ 10 11 G) accreting neutron stars, models for the disk-magnetosphere 
interaction based on this framework (see Ghosh & Lamb 1991 for a review) are generally consistent 
with the accretion torque behavior of the Be/X-ray pulsar transients, as observed by the Comp- 
fon/BATSE all-sky monitor (see Bildsten et al. 1997 and references therein). The Be/X-ray pulsar 
transients allow for a direct test of the predicted scaling of torque with accretion rate, since they 
span a wide range of accretion rates during their outbursts. However, the same models can only ac- 
count for the bimodal torque reversals observed in several persistent accreting pulsars if additional 
assumptions are introduced (Chakrabarty et al. 1997a, 1997b; Bildsten et al. 1997). Examples 
of suitable assumptions are a bimodal distribution of the mass transfer rate onto the pulsar, or 
a bimodal dependence on accretion rate of the orientation of the disk (Nelson et al. 1997; van 
Kerkwijk et al. 1998), of the orbital angular velocity of the accreting gas (Yi & Wheeler 1998), or 
even of the strength and orientation of any magnetic field produced in the disk (Torkelsson 1998). 

In weakly-magnetic (B <^ 10 10 G) accreting neutron stars, models of the disk-magnetosphere 
interaction have only been tested indirectly so far. Most neutron stars in low-mass X-ray binaries 
(LMXBs) show no periodic oscillations in their persistent emission. They have dipole magnetic fields 
^ 10 10 G, as inferred from their bursting (Lewin, van Paradijs, &: Taam 1995) and rapid variability 
behavior (see van der Klis 1998 for a review) as well as from their X-ray spectra (Psaltis k, Lamb 
1998). Their power-density spectra show various types of quasi-periodic oscillations (QPOs), and in 
particular the most luminous of the neutron-star LMXBs show ~ 15 — 60 Hz QPOs. These are called 
horizontal-branch oscillations (HBOs) and have frequencies that increase with mass accretion rate 
(van der Klis et al. 1985; van der Klis 1989). They have been interpreted as occurring at the beat 
frequency between the Keplerian frequency at the magnetic interaction radius (where the stellar 
magnetic field strongly interacts with the accretion disk) and the neutron star spin frequency (Alpar 
& Shaham 1985; Lamb et al. 1985). Models of the inner accretion disk and the disk-magnetosphere 
interaction can account for the observed scaling of HBO frequency with accretion rate if the neutron 
stars in these systems are near their magnetic spin equilibrium and if the inner accretion disks are 
radiation-pressure dominated (Psaltis et al. 1999b). 
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The discovery of ~ 200 — 1200 Hz QPOs (often occurring in pairs; hereafter kHz QPOs; see 
van der Klis 1998 and references therein) in the X-ray flux of many neutron-star LMXBs and the 
identification of the higher-frequency QPO with a Keplerian orbital frequency in the accretion disk 
(van der Klis et al. 1996; Strohmayer et al. 1996; Miller, Lamb, & Psaltis 1998) has introduced 
additional complications in studying the disk-magnetosphere interaction using HBO observations 
(see also Psaltis, Belloni, & van der Klis 1999a). The magnetic interaction radius inferred from 
the HBO frequencies is larger than the disk radius responsible for the higher-frequency kHz QPO. 
Therefore, if the higher-frequency kHz QPO is a Keplerian orbital frequency in the disk, then the 
magnetospheric beat-frequency model for the HBO can be valid only if a non-negligible fraction 
of the disk plasma is not threaded by the stellar magnetic field but remains in the disk plane 
inside the interaction radius (see Miller et al. [1998] and Psaltis et al. [1999b] for a discussion). 
This requirement has not yet been addressed in any theoretical model for the disk-magnetosphere 
interaction and therefore cannot be tested directly. 

The recent discovery of the first weakly-magnetic accretion-powered pulsar with the Rossi X- 
ray Timing Explorer (RXTE) allows the first direct test of disk-magnetosphere interaction models 
in the weak-field limit. SAX J1808. 4-3658 is a transient X-ray source that shows type I X-ray 
bursts (in't Zand et al. 1998) and coherent 401 Hz X-ray pulsations (Wijnands & van der Klis 
1998a), and is a member of a low-mass binary in a 2-hour orbit (Chakrabarty &; Morgan 1998). 
Its distance is ^4 kpc, as inferred from flat-topped type I X-ray bursts that are thought to be 
Eddington limited (in 't Zand et al. 1998), and its luminosity varied by ^2 orders of magnitude 
during the 1998 April/May outburst (Cui, Morgan, & Titarchuk 1998). 

In this paper we study the disk-magnetosphere interaction in SAX J1808. 4-3658. In §2, we com- 
pare the predictions of theoretical models for the inner accretion disk and the disk-magnetosphere 
interaction with observations of SAX J1808. 4-3658 and, in §3, we infer the magnetic field strength 
of the pulsar. In §4, we discuss our results and their implications for disk-magnetosphere interac- 
tion models and compare SAX J1808. 4-3658 to the millisecond radio pulsars and the non-pulsing 
neutron stars in LMXBs. 

2. DISK-MAGNETOSPHERE INTERACTION 

2.1. Assumptions and Formalism 

Throughout this paper we assume that most of the accreting gas around SAX J1808. 4-3658 
is confined in a geometrically thin accretion disk before interacting with the pulsar magnetic field. 
This assumption is supported by the transient nature of the source, which can be understood in 
terms of dwarf-nova-like accretion disk instabilities (see Chakrabarty & Morgan 1998). We also 
neglect the effect of wind mass loss from the inner accretion disk and of radiation drag forces, as 
well as all general relativistic effects. 
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The radius tq at which magnetic stresses remove the angular momentum of the disk flow and 
disrupt it can be estimated by balancing the magnetic and material stresses, (Ghosh & Lamb 1978, 
1979a) 

^4vrr 2 Ar = Mftr 2 , (1) 

where B p and B^ are the poloidal and toroidal components of the magnetic field, Aro is the radial 
width of the interaction region, M is the mass transfer rate through the inner disk, and Q(r) is the 
angular velocity of the gas at radius r. Assuming that the poloidal magnetic field is dipolar with 
magnetic moment fi and that the accretion flow is Keplerian, we obtain 

( 4 \ 1/7 

where 7b = (-B^/-B p )(Aro/Vo), M is the mass of the neutron star, and G is the gravitational 
constant. 

The boundary layer parameter 73 in equation (^) depends in general on all the other phys- 
ical quantities. Hence, different models for the inner accretion disk's environment and the disk- 
magnetosphere interaction generally predict different coefficients and scalings (Ghosh & Lamb 1992; 
see also Ghosh & Lamb 1978, 1979a; Wang 1996; Ghosh 1996). However, if the disk-magnetosphere 
interaction takes place in a region of the accretion disk where all physical quantities have a power- 
law dependence on radius (as is the case for most accretion disk models away from the stellar 
surface) and equation (0) describes angular momentum balance in the interaction region, then the 
Keplerian orbital frequency at ro can be written as 



^%,ohr 77^77—3 TV" , (3) 




where Mq is the solar mass; Me = 2 x 10 -8 M@ yr -1 is the Eddington critical accretion rate (at 
which the outward radiation forces in a spherically symmetric hydrogen flow balance gravity); and 
a, f3, 7, and z^k,o are parameters which depend upon the nature of the inner accretion disk. For 
our study, we used the parameters from four different inner disk models (see Table 1). We chose 
these models both because they span a wide range of physical conditions in the inner disk (one- 
and two-temperature plasmas, and gas and radiation-pressure-dominated flows) and a range of 
radiation processes, and because detailed calculations of the interaction with a stellar magnetic 
field have been performed for these models (see Table 1 and Ghosh & Lamb 1992 for a discussion 
of the physical processes involved in these models). 



2.2. Limits from the Detection of Coherent Pulsations 



We can use the presence of coherent X-ray pulsations in SAX J1808. 4-3658 over a wide range 
of mass accretion rate to test scaling (^) in the weak magnetic field regime. For a rotating neutron 
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star to appear as an accretion-powered pulsar, the stellar magnetic field must be strong enough to 
disrupt the Keplerian disk flow above the stellar surface^. Note that the orbital frequency uq at the 
interaction radius vq increases with mass accretion rate (a > 0). Therefore, the stellar magnetic 
field must be strong enough to disrupt the disk flow at radii larger than the neutron star radius 
R at the maximum mass accretion rate M max for which coherent pulsations were detected. This 
leads to a lower limit on the magnetic dipole moment, 

At the same time, the stellar magnetic field must be weak enough that accretion is not centrifugally 
inhibited at the minimum mass accretion rate M m i n for which coherent pulsations were detected. 
Therefore, the orbital frequency vq at this accretion rate must be larger than the spin-frequency of 
the neutron star v s . This leads to an upper limit on the magnetic dipole moment, 

\vK,o) \MqJ \ M E J V ; 

Of course, X-ray brightness modulations at the pulsar spin frequency may be possible even if 
one of the above two requirements is not met. For example, azimuthal variations in the efficiency 
of processes dependent on magnetic field strength (such as cyclotron emission, absorption, and 
resonant scattering) may be strong enough to produce X-ray pulsations, even if the magnetic field 
of the neutron star is dynamically unimportant and the accretion disk extends down to the stellar 
surface. However, if at any point during the outburst either the disk had reached the stellar 
surface or accretion had been centrifugally inhibited, then either the X-ray spectrum or the pulse 
amplitude should have changed abruptly. Instead, the X-ray spectrum of SAX J1808. 4-3658 was 
found to be remarkably stable (Gilfanov et al. 1998; Heindl & Smith 1998), and the pulse amplitude 
to increase only slightly from 4% to 7% as the inferred luminosity declined by more than two orders 
of magnitude during the outburst (Cui et al. 1998). 



Combining constraints (||) and (||) we obtain 

M > 0-047 ( (^T(J^) . 



M Q \ 401 Hz J \ M min J V 10 km 

Given the observed spin frequency and range of accretion rates at which coherent pulsations are 
detected, inequality (|6|) defines a limiting curve in the M-R parameter space for neutron stars. The 
scaling of this curve is cubic in R, independent of the details of the disk-magnetosphere interaction 



x Note that this is a necessary but not a sufficient condition: the X-ray brightness at infinity may not be modulated 
with any detectable amplitude if, e.g., the spin of the neutron star is perfectly aligned to its magnetic moment, or the 
neutron star is surrounded by a scattering medium that attenuates the pulsations (Brainerd & Lamb 1986; Kylafis 
& Phinney 1989). 
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Fig. 1. — Constraints on the 
neutron star mass and radius in 
SAXJ1808.4-3658. The solid lines 
represent mass-radius relations of 
stable neutron stars for various pro- 
posed equations of state (Cook et 
al. 1994). The dashed lines repre- 
sent the limits imposed by various 
models for the inner accretion disk 
(Table 1). The allowed region for 
each inner disk model is to the left 
of the corresponding dashed curve. 
The limiting curve for the 2S model 
essentially lies on the R = axis on 
this scale, and is not plotted. 
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model, with a coefficient that depends on a single model parameter^. As mentioned above, the 
X-ray spectrum of SAX J1808. 4-3658 was found to be remarkably stable during the outburst, so 
we can assume here that the ratio M max /M m i n scales as the ratio of the corresponding countrates in 
a given X-ray bandpass, independent of (unknown) bolometric corrections. For the 1998 April/May 
outburst of SAX J1808. 4-3658, this ratio was « 130 (Cui et al. 1998). 

Figure 1 compares the M-R relations for several proposed equations of state for stable neutron 
stars (see Cook, Shapiro, & Teukolsky 1994) to the limit (JgJ) for the four previously mentioned inner 
disk models applied to SAX J1808. 4-3658. Only the limiting curve for the 1R model, for which 
a ~ 0.2, is consistent with many proposed equations of state. In a thin accretion disk, radiation 
pressure dominates over gas pressure at radii (see Treves, Maraschi, & Abramowicz 1988) 

„, / r, \ /M\»/ R \-° M f M V' 76 
R« a5 ° 01 (55o)(s;) (iota) (uRlsj • (7) 

where r\ is a parameter in the range 10 2 -10 3 and a <^ 1 is the Shakura-Sunyaev viscosity parameter. 
During the decline phase of the outburst, the mass accretion rate was so small that r/R <C 1 
and the inner accretion disk in SAX J1808. 4-3658 was certainly not radiation-pressure dominated. 



2 While this work was in progress, we learned of a paper by Burderi & King (1998) in which a similar argument 
is used to constrain R by assuming a particular disk-magnetosphere interaction model (essentially the 1G model) 
and hence a specific value for a. However, the applicability of the 1G model to SAX J1808. 4-3658 is dubious, as we 
demonstrate here. In contrast, we use inequality (^) to constrain disk-magnetosphere interaction models and hence 
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Whether or not the disk was radiation-pressure dominated at the peak of the outburst is not clear 
given the uncertainty in the model parameter r\. However, it is very unlikely that there was a 
transition between a radiation-pressure and a gas-pressure-dominated inner accretion disk during 
the outburst, given the stability of the X-ray spectrum of the source and of the pulsed fraction of 
the emission (see discussion below). Radiation-pressure-dominated models for the inner accretion 
disk, like the one shown in Figure 1, are therefore not applicable to SAX J1808. 4-3658 since its 
inferred mass accretion rate is significantly sub-Eddington (in't Zand et al. 1998). 

From the three gas-pressure dominated models, the two-temperature ones (2B and 2S) are 
inconsistent with any of the proposed equations of state and are thus ruled out for SAX J1808.4- 
3658. The one-temperature gas-pressure dominated model (1G) is barely consistent with a very 
restricted range of M and R allowed by only a few equations of state. However, the maximum 
and minimum accretion rates at which RXTE detected coherent pulsations from SAXJ1808.4- 
3658 were determined by the peak luminosity of the outburst and the instrumental sensitivity, 
respectively, which are obviously unrelated to the limiting equalities in (Q) and (||). This makes it 
quite unlikely that the actual mass and radius of SAX J1808. 4-3658 lies very close to the limiting 
curve 1G in Figure 1. Therefore, either a remarkable coincidence has occurred or (more probably) 
the 1G model is also ruled out for SAX J1808. 4-3658. This is particularly interesting because the 
dependence of Vq on M predicted by the 1G model is very similar to the simple uq oc M 3 / 7 scaling 
law derived from equation (||) for constant 7b- This simple scaling law is widely used to describe 
the disk-magnetosphere interaction in accreting magnetic stars (and is in fact generally considered 
the standard model). Figure 1 implies that the scaling of v$ with M for SAX J1808. 4-3658 should 
be weaker than predicted by this simple scaling law. In particular, we find that a must be < 0.4, 
and hence 7b must be a function of M. 



2.3. Predicted Limits on the Accretion Torque in SAX J1808.4 3658 

An independent test of scaling (||) in the weak magnetic field regime is possible if the spin 
frequency derivative i> s in SAX J1808. 4-3658 due to accretion torques can be measured. For a 
Keplerian disk flow terminated at a radius ro, the accretion torque on the neutron star is given by 

2irI0 s = rjM (GMr ) l/2 , (8) 

where I is the neutron star moment of inertia and rj is a dimensionless quantity in which all 
the physics of the disk-magnetosphere interaction is parametrized (Ghosh & Lamb 1979b). For a 
prograde accretion disk, rj is positive and a strong function of both the magnetic field strength and 
the mass accretion rate (Ghosh & Lamb 1979b). For a retrograde accretion disk rj is negative and 
has an absolute value of order unity (Daumerie 1996). 

Applying constraints (0) and (||) to equation (||), we obtain an upper and a lower limit on the 
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spin frequency derivative for SAX J1808. 4-3658 predicted by accretion torque theory, 



1, ///10 45 gcm 
z> s < 1.8 x 10~ 13 ?? 1 



M/Mq J 
R s \ ( M \ - V3 



-1/3 



401 Hz, 



10 km/ \M Q 



0.06 M E 



(M min /M max ) Q / 3 
130V7 



and 



i q ///10 45 gem 2 ' 
z> s > 1.0 x lO" 13 /? - 7 



# s \ 3/2 / M Y l/2 ( M, 



0) 



(10) 



M/M Q y VlOkm/ VM y \0.06 M E/ 
The maximum observed 2-30 keV flux of SAX J1808. 4-3658 at the peak of the outburst was ~ 3 x 
10 9 erg cm -2 s _1 (Cui et al. 1998), which corresponds to a 2-30 keV luminosity of ~ 0.03(O e /47r)LE, 
where f2 e is the solid angle of emission and Le = 2 x 10 38 erg s _1 is the Eddington critical luminosity 
that corresponds to Me- However, the photon spectrum of SAX J1808. 4-3658 over most of the 
RXTE 2-100 keV bandpass is well described by a power-law of the form dN/dE oc E~ 2 , where E 
is the photon energy. The energy spectrum is thus very flat, vL v oc E°. As a result, the bolometric 
luminosity of the source must be larger than estimated above and is probably about twice the 
2-30 keV luminosity, given the expected upper- and lower-energy cut-offs of the spectrum (see also 
Heindl & Smith 1998 for spectral fits over the entire RXTE bandpass that show this effect). We 
therefore adopt M max = 0.06Me- 

Given a measurement of the spin frequency derivative in SAX J1808. 4-3658, we could use 
relations (||) and (|l0|) to place additional constraints on a for inner disk models or M and R for 
neutron stars, similar to the discussion above for the presence of coherent pulsations (Figure 1). 
However, variations of the spin frequency were not detected during the peak of the 1998 April/May 
outburst of SAX J1808. 4-3658, leading to an upper limit on the spin frequency derivative of \u s \ < 
7 x 10" 13 Hz s _1 (Chakrabarty & Morgan 1998). This measured upper limit is consistent with our 
predicted bounds on the spin frequency derivative @ and (|l(]) if rj ^ 7. 



3. THE MAGNETIC FIELD STRENGTH OF SAX J1808.4 3658 

In their discovery paper for SAX J1808. 4-3658, Wijnands & van der Klis (1998a) applied 
"standard" magnetospheric disk accretion theory (essentially the 1G model discussed in §2) and 
used the absence of centrifugal inhibition of accretion during the decline of the 1998 outburst to 
infer that B (2—6) x 10 s G. As the outburst declined, pulsations continued to be detected, leading 
to revised upper limits on the magnetic field strength [B <J (0.4 — 1.3) x 10 8 G, Cui et al. 1998; 
B £ few xlO 7 G, Gilfanov et al. 1998]. However, as we showed in §2.2, simple scaling arguments 
or models of the form (||) for the asymptotic regions of gas-pressure-dominated accretion disks 
cannot easily account for the coherent pulsations observed from SAX J1808. 4-3658 throughout its 
1998 outburst. Therefore, the previous estimates of the magnetic field strength of SAXJ1808.4- 
3658 which are based on similar scalings, are not valid. A more careful calculation is necessary 
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to infer the magnetic field strength in this system. Because no spin frequency derivative has been 
detected yet from this source (Chakrabarty & Morgan 1998), we cannot uniquely determine the 
stellar magnetic field strength, but can only constrain it. 

We begin with expression (^) for the interaction radius, which is derived directly from the 
angular momentum equation (|l]), and impose the same constraints that led to equations and 
(|5|). We set the maximum and minimum mass accretion rates at which coherent pulsations were 
detected to M max > 0.06Me and M m \ n < Af max /130 as inferred from observations (see Cui et 
al. 1998 and §2.3). We assume a neutron star radius in the 10-15 km range, as predicted by 
current equations of state (see Fig. 1). We also assume a neutron star mass in the 1.4-2.3 M & 
range, consistent with the inferred masses of recycled millisecond pulsars (Thorsett & Chakrabarty 
1999), the upper limits on neutron star masses in LMXBs inferred from kHz QPO observations 
(Miller et al. 1998; see also Zhang, Strohmayer, & Swank 1997), and the possible neutron star mass 
measurement in 4U 1820-30 (Zhang et al. 1998a). 

The limits on the magnetic field strength of SAX J1808. 4-3658 also depend upon the allowed 
range for the parameter 73, which in turn depends upon the fractional width Aro/Vo of the bound- 
ary layer and the magnetic pitch B^/Bp within this layer (see Ghosh & Lamb 1991 for a detailed 
discussion) . The fractional width should be significantly smaller than unity in order for the bound- 
ary layer approximation (|2|) to be valid. Assuming a dipolar poloidal field and neglecting mass 
loss from the disk as well as toroidal screening currents in both the disk and the magnetosphere, 
one finds Aro/ro — 0.3 (Ghosh &: Lamb 1991; Daumerie 1996). However, relaxing any of these 
assumptions leads to a significantly narrower boundary layer. We therefore conservatively adopt 
0.01 £ Ar/r < 1. 

The toroidal component Ba of the magnetic field is produced by the differential rotation 
of gas in the accretion disk with respect to the stellar spin. The magnetic energy stored in the 
magnetosphere increases with the twisting of the field lines (see Zylstra 1988 and references therein). 
There is an upper bound on this energy (and hence on the twisting of the field lines and the magnetic 
pitch) for a semi-infinite space in which the magnetic field strength goes to zero at infinity (Aly 
1984, 1991; Zylstra 1988). It has been argued that, above this upper bound, the only existing 
configurations are those with field lines that close at infinity (i.e., open field lines; Aly 1984, 1991). 
Simple estimates of the maximum twisting of the field lines (see Ghosh &: Lamb 1991 and reference 
therein) as well as detailed numerical calculations of the magnetospheric structure (Zylstra 1988) 
lead to maximum values of B^/Bp ~ 1. Combining this with our constraints on Aro/ro, we 
therefore adopt 0.01 ^ 7b(M) ^ 1- 

Applying all these constraints to SAX J1808. 4-3658, we obtain 

--^MM-(i) , ' , (^rfe) ,,i - ! (») 
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Fig. 2. — Constraints on the 
magnetic dipole moment fj, of 
SAXJ1808.4-3658 as a function of 
boundary layer parameter tb(M). 
The upper limit on fi comes from 
requiring that accretion is not cen- 
trifugally inhibited at the minimum 
M at which coherent pulsations 
were detected. The lower limit on 
fj, comes from requiring that the ac- 
cretion disk is terminated above the 
stellar surface at the maximum M 
at which coherent pulsations were 
detected. Note that 7b depends on 
M, and so must have different val- 
ues at the minimum and maximum 
M considered here. Conservatively, 
0.01 £ 7 B < 1. 
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which correspond to a dipolar magnetic field of a few times 10 8 G at the stellar pole. Figure 2 
shows the limits on the magnetic dipole moment as a function of the (unknown) parameter 7b (M) . 

Finally, we can estimate the mass accretion rate M eq for which the observed spin frequency of 
the neutron star is the equilibrium frequency (at which the accretion torque is zero) . The result is 



M, 



cq 



2 x 10 _11 o;r 7/3 



7B 
0.1 



A* 



10 26 Gem 3 / V1.4M 



M 



-5/3 



401 Hz 



7/3 



Mo yr- 1 , (13) 



where uj c is the critical fastness parameter (Ghosh Sz Lamb 1979b). The value of M eq depends 
sensitively on the (unknown) magnetic field strength and on tb(M). However, for values of these 
parameters consistent with the constraints imposed above, M eq agrees with the long-term average 
accretion rate of SAX J1808. 4-3658 inferred by Chakrabarty & Morgan (1998) from the fluence of 
the X-ray outburst and the outburst recurrence time. Moreover, in order for an LMXB with a 2 hr 
binary period to be below the separatrix between transient and persistent systems in the diagram 
of donor mass versus orbital period, it must have M <; 3 x 10~ n M Q yr -1 (van Paradijs 1996; 
King, Kolb, k, Szuszkewicz 1997), also consistent with the value estimated above (Chakrabarty & 
Morgan 1998). 
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4. 



DISCUSSION 



In this paper we have studied the disk-magnetosphere interaction in the first known accretion- 
powered, millisecond pulsar, SAX J1808. 4-3658. We have demonstrated that various simple models 
of inner accretion disks are not consistent with observations. We have also argued that the magnetic 
field strength of SAX J1808. 4-3658 is a few times 10 8 G at the stellar pole, using very general 
constraints on the properties of the neutron star and on the physics of the inner accretion disk flow. 
In this section, we discuss the implications of our results for models of the disk-magnetosphere 
interaction around weakly magnetic neutron stars. We also compare SAX J1808. 4-3658 to the 
recycled millisecond radio pulsars and to the non-pulsing neutron-star LMXBs. 

4.1. The Disk-Magnetosphere Interaction in Weakly Magnetic Neutron Stars 

BATSE observations of strongly magnetic (B ~ 10 12 G) accretion-powered pulsars in Be/X-ray 
transients have shown that the scaling of the radius ro of interaction between the stellar magnetic 
field and a gas-pressure dominated disk flow can account for the observations (Ghosh 1996; Finger, 
Wilson, & Harmon 1996; Bildsten et al. 1997). However, in §2 we showed that this same scaling 
is inconsistent with the detection of coherent pulsations throughout the 1998 April/May outburst 
of SAX J1808. 4-3658, for most equations of state of neutron star matter. This is not surprising, 
given the very different physical conditions in the interaction regions around weakly- and strongly- 
magnetic neutron stars. 

In writing the stress balance equation ([l]) and the scaling we made a number of implicit 
assumptions regarding the properties of the neutron star and the inner accretion disk flow: first, 
that the stellar magnetic field is dipolar; second, that the pulsar magnetic moment is parallel to its 
spin axis; third, that the dominant mechanism for removing angular momentum from the accreting 
gas in the disk is magnetic stress; and finally, that the gravitational field is Newtonian everywhere 
and hence that stable, circular Keplerian orbits exist at all radii. Although these assumptions are 
valid (or at least are reasonable approximations) in the case of a strongly magnetic neutron star, 
they break down when ro is comparable to the stellar radius (see Lai 1998 for solutions of the 
structure of the inner accretion disk where some of these effects have been taken into account in a 
simple way). 

Let us examine the validity of each of these assumptions. First, magnetic stresses would remove 
angular momentum faster than predicted by equation (^) if higher-order multipoles were present 
in the pulsar magnetic field. For example, if the stellar magnetic field could be described entirely 
by a multipole of order I (i.e, the magnetic field strength in the equator was B = /^/r' +1 ), then 
the Keplerian frequency at ro would be 




(14) 



Equation (|l4|) shows that, if 7b depends only weakly on M, the dependence of vq on M weakens 
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Fig. 3. — The effect of the presence of a 
quadmpole magnetic moments on the con- 
straints on the neutron star mass and ra- 
dius in SAXJ1808.4-3658. The solid lines 
represent mass-radius relations of stable 
neutron stars for various proposed equa- 
tions of state (Cook et al. 1994). The 
dashed lines represent the limits imposed 
by the 1G model for the inner accretion 
disk and for different fractional contribu- 
tions of an aligned quadrupole moment to 
the magnetic field of the neutron star. The 
allowed region for each model is to the left 
of the corresponding dashed curve. 
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as the order I of the multipole increases: for a dipolar (I = 2) field vq oc M 0A3 , for a quadrupolar 
(I = 4) field vq oc M ' 2 , etc. As a result, the existence of intrinsic or induced multipole moments in 
the magnetosphere can significantly weaken the M-dependence of vq, hence making the scaling ( |l4|) 
consistent with observations of SAX J1808. 4-3658. As an illustration of this, Figure 3 shows the 
constraints imposed by the 1G model on the mass and radius of the neutron star in SAX J1808.4- 
3658 for various contributions of an aligned quadrupole moment to the pulsar magnetic field. 
Clearly, even a modest quadrupole moment (which would be completely negligible at distances 
larger than a few stellar radii from the surface) is enough to make the 1G model consistent with 
observations. 

The intrinsic multipole moments in neutron star LMXBs are thought to be weaker than their 
dipole moments, based on the current spin frequencies and spin frequency derivatives of the millisec- 
ond radio pulsars thought to be their descendants (Arons 1993). However, electrical currents in the 
accretion disk and the magnetosphere can significantly enhance the strength of multipole moments 
by altering the magnetic field topology in the interaction region (see, e.g., Ghosh & Lamb 1979a). 
Moreover, electrical currents on the neutron star induced by the presence of an accretion disk can 
produce multipoles as strong as the intrinsic dipole, provided that the disk extends very close to 
the stellar surface (Psaltis, Lamb, & Zylstra 1996). The corotation radius r co = (GM/Att 2 ^) 1 ^ 3 
(which is an upper bound on tq when coherent pulsations are detected; see, e.g., Ghosh & Lamb 
1979a) is only ~ 2.8 times the neutron star radius in SAX J1808. 4-3658, and therefore none of the 
above effects are negligible. 

Second, if any of the magnetic moments of SAX J1808. 4-3658 is not aligned to the stellar 
rotation axis (which is actually a necessary condition for the system to appear as a pulsar) then 
the scaling of the inner Keplerian disk radius with accretion rate may be different than what is 
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Fig. 4. — Inferred pulsar dipole mag- 
netic field strength at the stellar equator, 
as a function of spin period. The small 
dots are normal radio pulsars (Taylor et 
al. 1993). The open circles are millisec- 
ond radio pulsars which are not members 
of globular clusters, and for which the in- 
ferred field strengths have been corrected 
for the apparent period derivative caused 
by their transverse velocities (Camilo et al. 
1994). The error bar shows the allowed 
range of field strengths for SAXJ1808.4- 
3658. The spin-up line (solid) is shown 
here for illustrative reasons only and rep- 
resents the minimum period of recycled 
pulsars for a given field strength, assum- 
ing a specific model for the inner accretion 
disk and the disk-magnetosphere interac- 
tion. The death line (dashed) is an em- 
pirical estimate of the maximum period at 
which pulsars have detectable radio emis- 
sion for a given field strength. 

predicted by the models considered above. This case has never been treated in the literature 
because it corresponds to a time-dependent non-steady-state problem. It is beyond the scope of 
the current paper to generalize models of the interaction between the pulsar magnetic field and the 
accretion disk to time-dependent situations. We will therefore not consider this case any further. 

Third, when the accretion disk penetrates very close to the neutron star, the effect of magnetic 
stresses is amplified by radiation drag forces that can efficiently remove angular momentum from 
the accreting gas (Miller & Lamb 1996; Miller et al. 1998). Finally, the specific angular momentum 
of gas at radii comparable to the innermost stable orbit around a compact object has a weaker 
dependence on radius than estimated in a Newtonian gravitational field. 

All the effects discussed above result in a weaker M-dependence of vo, thus making scaling (|2|) 
consistent with observations of SAX J1808. 4-3658. 

4.2. Comparison with Millisecond Radio Pulsars 

When mass transfer onto the neutron star ends, SAX J1808. 4-3658 may appear as a millisecond 
radio pulsar (Wijnands & van der Klis 1998a; see also Bhattacharya & van den Heuvel 1991). 
Figure 3 compares the current spin frequency and inferred magnetic field strength of SAX J1808.4- 
3658 at the stellar equator with those of the known normal and millisecond radio pulsars (Taylor, 
Manchester, & Lyne 1993). We only include millisecond radio pulsars in the Galactic disk, since the 
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inferred spin frequency derivatives (and hence the magnetic field strengths) of pulsars in globular 
clusters are probably contaminated by gravitational acceleration in the potential of the cluster. 
Moreover, we use magnetic field strengths that have been corrected for the effect of the apparent 
spin frequency derivative caused by the non-negligible transverse velocities of the pulsars (Camilo, 
Thorsett, & Kulkarni 1994; Toscano et al. 1999). Note, however, that the inferred magnetic dipole 
moments assume orthogonal dipole rotators, and model-dependent systematic effects have not been 
taken into account. 

The mass of the companion of SAX J1808. 4-3658 is <^ O.3M , and hence the binary must be 
close to the termination of its X-ray phase (Chakrabarty & Morgan 1998). The current magnetic 
field strength and spin frequency of the neutron star should thus be very similar to that of the 
descendant millisecond radio pulsar. Figure 3 shows that the descendant of SAX J1808. 4-3658 
should appear as a normal millisecond radio pulsar in the Galactic disk. Given the small number 
statistics and the uncertainty in our estimate of the magnetic field strength in SAX J1808. 4-3658 
and in the millisecond radio pulsars, it is not possible to determine whether the field strength 
of SAX J1808. 4-3658 is lower or higher than those in the known millisecond radio pulsars. Such 
a determination would be of interest in comparing SAX J1808. 4-3658 to the other non-pulsing 
LMXBs. 



4.3. Comparison with Non-Pulsing LMXBs 

SAX J1808. 4-3658 is the only known weakly-magnetic accreting neutron star with persistent 
coherent pulsations at the neutron star spin frequency. Upper limits on the amplitudes of periodic 
oscillations in other neutron-star LMXBs range from <; 1% for the most luminous Z sources to a 
few percent for the less luminous atoll sources (see Vaughan et al. 1994 and references therein). 
Coherent millisecond oscillations, however, possibly at the spin frequencies of the neutron stars, 
have been observed in several neutron-star LMXBs during Type I X-ray bursts (Strohmayer, Swank, 
k Zhang 1998). 

The absence of persistent coherent pulsations in most LMXBs might be due to neutron star 
magnetic field strengths that are significantly different than that in SAX J1808. 4-3658. Specifically, 
the field strengths in the non-pulsing LMXBs may be so strong that accretion is centrifugally 
inhibited, or so weak that the accretion disk extends all the way to the neutron star surface. 
Alternatively, the neutron stars in the non-pulsing LMXBs may be surrounded by optically thick 
scattering media that spread the pulsar beams and attenuate the oscillations below present detection 
thresholds (Brainerd & Lamb 1987; Kylafis & Phinney 1989). We consider here each possibility 
separately. 

It is rather unlikely that the field strengths of all non-pulsing LMXBs are so strong that accre- 
tion is centrifugally inhibited, since at least some are accreting at near-Eddington rates. Moreover, 
if the field strengths were all so strong, then the observational fact that all known millisecond radio 
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pulsars lie below the so-called spin-up line in the B-P diagram (Fig. 3) would have to be a selection 
effect or a coincidence. This is because the spin-up line defines the minimum spin period of a 
recycled pulsar (which is the equilibrium period at the Eddington accretion rate) for a given field 
strength, provided that accretion is not centrifugally inhibited. 

On the other hand, if the field strengths in the non-pulsing LMXBs are significantly weaker 
than that of SAX J1808. 4-3658, then a significant fraction of the gas in the accretion disk may 
be reaching the stellar surface in the disk plane, producing undetectable periodic oscillations at 
the neutron star spin frequency. This is apparently consistent with identifying the frequencies of 
the observed kHz QPOs in the non-pulsing LMXBs with Keplerian orbital frequencies in the inner 
accretion disk (van der Klis et al. 1996; Strohmayer et al. 1996; Miller et al. 1998). 

However, there are at least two other transient LMXBs (Aql X-l and 4U 1608—52) whose 
outburst histories, X-ray spectral properties, and broad-band noise features are similar to those of 
SAX J1808. 4-3658 (see, e.g., Heindl & Smith 1998; Wijnands & van der Klis 1998b), and yet do 
not show periodic pulsations in their persistent emission. At the same time, these systems do show 
kHz QPOs and have other timing and spectral properties characteristic of the non-pulsing LMXBs 
(Zhang, Yu, & Zhang 1998b; Yu et al. 1997). In particular, the X-ray spectrum of 4U 1608-52, 
which is a direct probe of the structure and dynamics of the inner accretion disk flow, is nearly 
identical to that of SAX J1808. 4-3658 (Heindl &: Smith 1998), with the same power-law index and 
high-energy cut-off. It seems unlikely that the X-ray spectra of these two sources could be so similar 
if their field strengths (and thus their inner disk radii) differ significantly. This is particularly true 
of the spectra at high energies (20-100 keV), which must be formed close to the neutron star surface 
and are thus highly sensitive to the inner disk flow geometry. In fact, if the high-energy spectra are 
produced by Comptonization of high-harmonic self-absorbed cyclotron photons (Psaltis et al. 1995), 
then the stellar field strengths can be inferred from the hardness of the 2 — 20 keV X-ray spectra 
(Psaltis 1998; Psaltis & Lamb 1998). Indeed, the magnetic field strengths of the low-luminosity 
atoll sources, such as Aql X-l and 4U 1608—52, have been inferred to be ~ 5 x 10 s G (see Psaltis 
& Lamb 1998), i.e., very similar to the magnetic field strength of SAX J1808— 36 we estimated in 
§3- 

If the field strength in SAX J1808. 4-3658 is comparable to that of some non-pulsing LMXBs, 
then the absence of periodic oscillations in the persistent emission from most LMXBs must be 
attributed to another, external effect. The presence of a scattering medium around the neutron 
star and its magnetosphere has been invoked previously as an explanation of the non-detection of 
periodic oscillations in neutron-star LMXBs (Brainerd & Lamb 1987; Kylafis & Phinney 1989). The 
high (r ~ 4 — 5) optical depths inferred from X-ray spectra of the Z and atoll sources are consistent 
with this hypothesis (see Psaltis et al. 1995 and references therein). However, RXTE /HEXTE 
observations of SAX J1808. 4-3658 indicate an optical depth of r = 4.0 ± 0.2, based on the slope 
and high-energy cutoff of the X-ray spectrum, at the same time that pulsations were detected 
(Heindl & Smith 1998; Heindl 1998, private communication; note, however, that Gilfanov et al. 
[1998] analyzed the same data and reached a different conclusion.). This is very similar to the 
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optical depth inferred for the non-pulsing source 4U 1608—52 (Heindl & Smith 1998), which we 
argued above has a comparable magnetic field strength. Thus, if the optical depth of SAX J1808.4- 
3658 is indeed ;> 3, then the reason coherent pulsations are detected only from this LMXB cannot 
simply be weaker attenuation of the beaming oscillation in the scattering medium. 

A number of additional effects may be responsible for the detectability of coherent pulsations 
in SAX J1808. 4-3658. For example, its very small measured mass function is suggestive of a nearly 
face-on viewing angle of the binary (Chakrabarty & Morgan 1998). In this case, the pulsed emission 
may propagate through a region of lower optical depth close to the pulsar rotation poles, whereas 
the majority of the (non-pulsed) X-ray photons may travel through optically thick regions nearby. 
However, it is important to note that the presence of a shallow partial eclipse in SAX J1808. 4-3658 
has been suggested (Chakrabarty & Morgan 1998; Heindl & Smith 1998). Confirmation of this 
feature would rule out a face-on orientation of the binary, thus challenging the above explanation. 
Alternatively, different stellar field topologies in SAX J1808. 4-3658 and the non-pulsing neutron 
star LMXBs, possibly related to different prior evolution or even to different magnetic inclinations, 
may be responsible. 

A number of other LMXBs share many characteristics of SAX J1808. 4-3658 and are thus 
good candidates for detecting periodic pulsations in their persistent emission. For example, as 
we discussed above, Aql X-l and 4U 1608—52 are transients with X-ray spectra very similar to 
SAX J1808. 4-3658, suggesting similar field strengths and inner disk flows. The nearly sinusoidal 
optical orbital-phase lightcurves of three other LMXBs (4U 1636-53, GX 9+9, and 4U 1735-44) 
suggest that these are also viewed nearly face-on (see, van Paradijs & McClintock 1995). Finally, 
GS 1826-26 (Homer, Charles, & O'Donoghue 1998) and MS 1603+2600 (Hakala et al. 1998) arc 
other LMXBs with low-mass companions in ~ 2-hour orbits, in which mass transfer is probably 
driven by angular momentum losses due to gravitational radiation. These binaries are thus in 
an evolutionary stage very similar to that of SAX J1808. 4-3658. In particular, GS 1826—26 is 
a transient X-ray burster with a low- amplitude modulation of its optical brightness, again possi- 
bly suggesting a face-on orientation (Homer et al. 1998). Detecting or imposing stringent upper 
limits on coherent pulsations in the persistent emission of any of these sources will be crucial in 
understanding the characteristics of SAX J1808. 4-3658 and the relationship between LMXBs and 
millisecond radio pulsars. 
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Table 1. Scaling of the Keplerian Frequency at the Coupling Radius a 



Disk Modcl b 


VKfl (Hz) 


a 


a 


7 


1G 


430 


0.38 


-0.87 


0.82 


1R 


210 


0.23 


-0.77 


0.70 


2B 


80 


0.72 


-0.86 


0.43 


2S 


50 


2.55 


-1.20 


-0.60 



a After Ghosh & Lamb 1992. 

b lG: Optically thick, gas-pressure dominated (GPD) disk; 1R: Optically thick, 
radiation-pressure dominated (RPD) disk; 2B: Two-temperature, optically thin 
GPD disk with Comptonized bremsstrahlung; 2S: Two-temperature, optically 
thin GPD disk with Comptonized soft photons (for references see Ghosh & Lamb 
1992). 



